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Abstract 

A one-dimensional tight-binding model of electrons with long-range Coulomb in- 
teractions is studied in the limit where double site occupancy is forbidden and the 
Coulomb coupling strength V is large with respect to the hopping amplitude t. The 
quantum problem of a kink-antikink pair generated in the Wigner lattice (the clas- 
sical ground state for t = 0) is solved for fillings n = where s is an integer larger 
than 1. The pair energy becomes negative for a relatively high value olV , Vc/t ~ . 
This signals the initial stage of the quantum melting of the Wigner lattice. 
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1 Introduction 



Three decades ago, Michael Rice and collaborators have introduced nonlin- 
ear phase excitations of a charge-density-wave condensate as a new type of 
charged states and referred to them as ip particles [1]. Their arguments were 
based on the Peierls instability of a one-dimensional coupled electron-phonon 
system, in the limit where the charge-density wave amplitude is much smaller 
than the average electronic charge density. Shortly after, Hubbard discussed 
qualitatively the role of dimer pairs (domain boundaries) in a nearly quarter- 
filled band with dominant long-range Coulomb interactions [2]. As we will 
argue, the two types of charged states are closely related. 

We consider a one-dimensional system of spinless fermions described by the 
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Hamiltonian 

H ^-t ^{cjci+i + 4+iQ) + ^ ^ ViTliTli+i , (1) 
i i 1>1 

where cj and q are fermionic creation and annihilation operators, respectively, 
the occupation number of site i and Vi = V/l represents the 
long-range Coulomb potential. Such a model may describe molecular chain 
compounds where the interaction between two (valence) electrons on the same 
molecule is so large that double occupancy can safely be discarded. The spin 
quantum number is then redundant because the exchange of two electrons is 
dynamically forbidden. 

We have in mind systems where the average site occupation is a rational 
number between and 1, n — r/s. Depending on the band filling the fermions 
then represent either electrons or holes. The classical ground state (t = 0) is 
charge-ordered and forms a "generalized Wigner lattice" [2] . In the special case 
n = 1/s the unit cell contains s sites, one of which is occupied. A generalized 
Wigner lattice is clearly insulating and therefore the long-range Coulomb 
interactions remain unscreened. Insulating charge-ordered phases have been 
reported for organic chain compounds [3] and in several cases the nature of 
charge ordering is more likely that of a Wigner crystal than that of a (small 
amplitude) charge-density wave. 

We have recently used a variational wave function for describing the effects 
of a small but finite hopping term (t <^ V) [4]. Here we consider the same 
regime, but follow a different route by examining the stability of the classical 
ground state (the generalized Wigner lattice) with respect to the spontaneous 
generation of charge defects, as was done in ref. [5] in the case of electron- 
phonon interactions. Such defects occur in pairs (kinks and antikinks) and 
can be viewed as fractionally charged particles, in close analogy to the (p 
particles of Michael Rice. 



2 Kink— antikink pair for n — 1/2 

We start our discussion for a density n = 1/2. In order to guarantee overall 
charge neutrality we introduce a rigid compensating background (n = 1/2 
charge of opposite sign at each site). The classical ground state in this case 
corresponds to alternating filled and empty sites (Fig. la). Note that this con- 
figuration, already at the classical level, is doubly degenerate, since exchanging 
the empty and occupied sites does not change the overall energy. The most 
simple defect in this perfect structure, called kink in the following, is a domain 
boundary separating the two possible ground state configurations. It can be 
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a) x»x»x»x»x» classical ground state 

b) x»»xx»x»x« kink-antikink pair at d=l 

c) x»»x«xx«x» kink-antikink pair at d=2 

d) X • X • x • • kink-antikink pair at d=3 



Fig. 1. Classical ground state and low-energy defect pairs at n = 1/2. Dots/crosses 
stand for filled/empty sites respectively. The kinks and anti-kinks correspond to 
the shaded areas. The arrows indicate the hopping processes leading to a given 
configuration. 

visualized as a pair of occupied nearest neighbor sites (a dimer) and is nec- 
essarily accompanied by an antikink, a pair of empty nearest-neighbor sites. 



The classical energy of a kink-antikink pair is readily evaluated. Introducing 
a pair of size 2d, i.e. a kink and an antikink separated by a distance 2d (see 
Fig. 1), costs an energy 

d oo 

AW = E E [^2m+l + V^m-l - 2V2m]. (2) 
p=l m=p 



This sum can be transformed into a power series in l/{2d) with leading terms 
V V 1 VI 



Here the first term represents the creation energy of two well separated defects. 
The second term is the Coulomb attraction between kink and antikink with 
a coupling reduced by a factor of 4 as compared to the original Hamiltonian 
(1). Therefore the defects carry effective charges ±1/2 (measured in units of 
the electronic charge). This result agrees with what one obtains by adding or 
subtracting an electron to the classical ground state configuration [2]. In fact, 
adding an electron at an empty site results in a high-energy configuration 
with a three-electron cluster. The energy is lowered by moving away one of 
these electrons, thus creating a pair of defects (two kinks in our language). 
Since the overall added charge is -1, each kink carries a charge -1/2. Similarly, 
removing an electron produces two antikinks, each of which carries a charge 
+1/2. The third term in Eq. 3 can be interpreted as the interaction of two 
dipole moments pointing in opposite directions parallel to the chain axis. The 
size of the dipoles is equal to the fractional charge times half a lattice constant, 
in our units 1/4. We conclude that kinks and antikinks behave like fractionally 
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charged particles with electric dipole moments. The dominant interaction is 
the Coulomb attraction, even at the shortest possible distance 2d — 2. 

We now turn to the quantum problem of a defect pair, i. e. we diagonalize the 
Hamiltonian within the subspace of states \m,d), corresponding to a kink at 
site m — d and an antikink at site m + d {m thus indicates the center-of-mass of 
the pair). In general the effect of a hopping event on the state |m, d) is either 
to move an existing defect by two lattice sites, which modifies both the size 
{2d — > 2d ± 2) and the center {m ^ m — ±1) of the pair (see Fig. 1), or to 
create (annihilate) an additional pair. Restricting ourselves to the subspace of 
single-pair states, creation and annihilation processes are forbidden and the 
quantum problem reduces to the eigenvalue equation 

[H - A(c?)] \m,d) + t [\m - l,d + I) + \m + I, d - 1) 

+ |m + l,(i+l) + |m-l,d-l)] = . (4) 

The center-of-mass and relative motions can be separated by introducing the 
Bloch superposition 

IV') = E e'^^m\m.d) . (5) 

m,d>l 



Introducing this Ansatz into Eq. (4) leads to the following eigenvalue equation 
for the wave function ip{d), 

Eij{d) = A{dyip{d)-2t cos K[tp{d+ I) +'tp{d-l)] , d>l. (6) 



Seeking for the lowest-energy wave function we restrict ourselves to the pair 
at rest {K — 0). The classical ground state (no defect pairs) is excluded 
by imposing the boundary condition "0(0) = 0. This eigenvalue problem can 
be treated to arbitrary accuracy numerically, but more insight is gained by 
solving it in the continuum limit, tjj{d,) 'ip{x), ip{d± 1) iIj{x±2), where x 
is the position in units of the distance between neighboring sites. Neglecting 
the dipolar interaction in Eq. (3) we obtain the eigenvalue equation 



where e — E — V/2 + At. Eq. (7) is equal to the radial eigenvalue equation 
of the hydrogen problem in the s-wave channel, and the boundary conditions 
are also the same, ip{0) = 0. Therefore we can immediately write down the 
lowest-energy wave function of a pair, 

V'o(a;) = 2ao^^^a;e-^/"°, x>0 , (8) 
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with an effective Bofir radius 

00 = -^ (9) 



and an energy 

V 



When the interaction strength is lower than the critical value Vc/t — 8.27, 
Eq < and the formation of a kink-antikink pair becomes favorable by virtue 
of the gain in energy achieved through delocalization as a Bloch state. ^ 

At the critical interaction strength the Wigner lattice will be destabihzed due 
to a proliferation of kink-antikink pairs, but we do not expect a true transition 
to a one-dimensional metallic state (a Luttinger liquid), because for arbitrarily 
small V the system is unstable with respect to a charge-density wave with the 
same period as for the Wigner lattice [6]. Therefore the critical strength Vc 
indicates a crossover between a Wigner lattice and small-amplitude charge- 
density wave. Note that in the crossover region the binding energy V^/512t 
makes a negligible correction to the quantum mechanical energy of a pair: 
the main contributions to Eq. (10) are the delocalization energy — 4t, and the 
formation energy of the isolated defects V/2. 



3 Kink-antikink pairs for n — 1/s 



We extend now the considerations of the previous section to other filling fac- 
tors n 7^ 1/2. For a general value of n already the classical limit {t — 0) is quite 
involved, leading to complicated charge patterns [2]. Therefore it is worthwhile 

to make contact to the opposite limit, V < t, where a simple picture is avail- 
able. Mean-field theory apphed to the Hamiltonian (1) yields a ground state 
with a modulated charge density 

{5ni) = a cos{Ql + (pi) , (11) 



where Q = 2kF = 27rn in the case of spinless fermions, and a <^ n for V < t. 
The phase is an arbitrary constant for an incommensurate situation {n 

^ The continuum description is appropriate provided that the extent (,t) = 96t/V is 
larger than the average distance between particles, a condition that is well satisfied 
at the critical coupling Vc- Indeed, solving the original discrete problem of Eq. (6) 
numerically yields Vc/t = 8.25, which is very close to the continuum estimate. 
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a) "xxx^xxx^xxx^xxx* classical ground state 

b) "xxx^xx^xxxx^xxx* pairatd=l 

c) -xxx-xx-xxx^xxxx* pairatd=2 

Fig. 2. Classical ground state and low energy defect pairs at n = 1/4. The symbols 
are the same as in figure 1. 

irrational), but locked to one of s possible values in the commensurate case, 
where n = r/s with integer numbers r and s [7]. Adding an electron to an 
incommensurate charge-density wave results in a slight shift of the wave vector 
Q, Q — i> [1 + {1/L)]Q, where L is the chain length. In view of Eq. (11), this 
effect can also be attributed to a phase shift (fi = 2711/ L. For a commensurate 
case the phase achieves the change of 2tt from / = to / = L through s steps, 
i.e. there are s y^-particles, each of which carries a fractional charge — 1/s. 

Let us now return to the large V limit and restrict ourselves to simple ratios 
n = 1/s, s > 3. In this case the classical ground state is s-fold degenerate and 
there are several types of domain walls separating the different configurations. 
It turns out that the low-energy domain walls are those which connect nearby 
ground state configurations, i.e. Wigner lattices where the locations of elec- 
trons differ by one lattice constant (see Fig. 2). These "kink" and "antikink" 
defects again occur in pairs. Repeating the arguments of the steps of the pre- 
vious section, we evaluate the classical energy of a kink-antikink pair of size 

d oo 

^{d) = J2T. (Vsm+l + Km-l - 2Vsm) ■ (12) 
p=l m=p 



The summation can again be performed for large sizes, and we find up to first 
order in 1/ (sd) 



V 

AW = - 



TT 



1 cot - 



TT 



V 1 

sd 



(13) 



The first term, the pair creation energy A^o at infinite separation, agrees with 
the corresponding quantity in Eq. (3) for s = 2. The second term corresponds 
to the Coulomb attraction for two particles with fractional charges ±l/s at 
a distance sd, in agreement with the counting argument presented above for 
the v?-particles. 

We can again use the continuum limit for calculating the lowest -energy quan- 
tum pair state for general s. The wave function is again given by Eq. (8), but 
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with a Bohr radius 

«o = ■ (14) 



Thus the pair size increases strongly with decreasing density n — 1/s. Corre- 
spondingly, the binding energy decreases, as seen in the pair energy 

Eo^-At + A^-^^, (15) 



where the first two terms dominate for large s. In fact, for s ^ 1 we can safely 
use the asymptotic value 

Eo^-U + ^V (16) 



to estimate the instability point, where this energy becomes negative. The 
result, 

(17) 

T^^^ ^^^^ 



agrees well with our previous variational estimate [4]. This is not surprising 

because the creation of a pair can be achieved through a hopping event that 
moves a particle out of the ground state configuration. The probability of 
such a hopping event was evaluated in Ref. [4], and was used to determine a 
criterion for the instability of the generalized Wigner lattice. 

In contrast to the case n = 1/2, for lower density this instability does not signal 
a crossover to a small-amplitude charge-density wave, it rather indicates that 
part of the electronic charge is spilled over to neighboring sites of the classical 
Wigner lattice. 



4 Discussion 



In this paper we have determined the lowest-energy quantum state of a kink- 
antikink pair in a one-dimensional generalized Wigner lattice. We have cal- 
culated the critical value Vc of the interaction strength below which the pair 
energy is negative, and charge defects will be generated spontaneously. Kinks 
(or antikinks) are the strong-coupling analogs of the (^-particles studied a 
long time ago by Michael Rice and collaborators [1]. 
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The true quantum ground state of electrons interacting through long-range 
Coulomb forces contains kink-antikink pairs due to quantum fluctuations, 
even above Vc- It would be interesting to proceed from the single-pair solu- 
tion to that of an arbitrary number of pairs. This step is highly non-trivial, 
although at first sight it looks similar to that from the Cooper problem to the 
BCS wave function. One of the difficulties arises from the unknown statistics 
of kinks and antikinks, another from their non-local character. 

If the density of pairs in the ground state is large enough, they loose their 
identity, and kinks and antikinks may move rather independently. In this case 
one can imagine a d.c. charge transport due to moving defects. Whether such 
a mechanism is responsible for the observed Drude peak in the Bechgaard 
salts [8] is an interesting open question. 
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